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^Vbstract 



i— We present theoretical work in which the degree of electrostatic 
t^^ oupling across a charged lipid bilayer in aqueous solution 
Os analyzed on the basis of nonlinear Poisson-Boltzmann the- 
^fry. In particular, we consider the electrostatic interaction of a 
jingle, large macroion with the two apposed leaflets of an op- 
positely charged lipid bilayer where the macroion is allowed 
io optimize its distance to the membrane. Three regimes are 
i-^dentified: a weak and a high macroion charge regime, sep- 
arated by a regime of close macroion-membrane contact for 
intermediate charge densities. The corresponding free energies 
are used to estimate the degree of electrostatic coupling in a 
lamellar cationic lipid-DNA complex. That is, we calculate to 
^Hvhat extent the one- dimensional DN A arrays in a sandwich-like 
rfipoplex interact across the cationic membranes. We find that, 
Vln spite of the low dielectric constant inside a lipid membranes, 
Zlh ere can be a significant electrostatic contribution to the ex- 
^jjerimentally observed cross-bilayer orientational ordering of 
f~tfie DNA arrays. Our approximate analytical model is comple- 
\Qnented and supported by numerical calculations of the electro- 
<^^tatic potentials and free energies of the lamellar lipoplex ge- 
~~^%metry. To this end, we solve the nonlinear Poisson-Boltzmann 
Equation within a unit cell of the lamellar lipoplex using a new 
lattice Boltzmann method. 

introduction 

•^ipids self-assemble in aqueous solution into a variety of mor- 
. £j>hologically different structures among which the planar bi- 
^jayer is the biologically most relevant one. Almost all natu- 
jjfally occurring lipid membranes contain to some degree acidic 
lipids which forms the basis for electrostatically-driven protein 
adsorption. In fact, the adsorption of macroions onto oppositely 
charged lipid membranes is an ubiquitous phenomenon not 
only in the living cell but is also common in technological and 
pharmaceutical applications (Evans and Wennerstrom, 1994). 
For example, the formation of cationic lipid-DNA complexes 
( |de Lima et ah, 2001) - one of the most promising carriers for 
non-viral gene delivery - is based on electrostatic interactions 
and subsequent complexation of cationic lipids and (the nega- 
tively charged) DNA. 

What determines the electrostatic properties of a lipid bi- 
layer is not only the mole fractions (and valencies) of the 
charged lipid species but also the low dielectric constant of the 
membrane interior ( Israelachvili, 1992 1. The large drop in the 



dielectric constant across a lipid bilayer makes it energetically 
unfavorable for large electric fields to form and thus dramat- 
ically diminishes the electrostatic coupling between the two 
charged monolayers of a lipid bilayer. Similarly, the electro- 
static interaction between two charged macroions, adsorbed on 
the opposite leaflets of a lipid bilayer, is masked. 

We mention two systems where the electrostatic inter- 
action across a lipid bilayer may be significantly enhanced: 
The first concerns the deposition near a lipid membrane of 
an anorganic crystal, calcium pyrophosphate dihydrate, lead- 
ing to crystal-induced membranolysis. A recent Molecu- 
lar Dynamics simulation of crystal deposition on the zwit- 
terionic membrane dimyristoylphosphatiadylcholine (DMPC) 
suggests that membrane rupture starts in the apposed mono- 
layer and is induced by electrostatic interactions across the 
bilayer (Dalai et ai, 2006 1. The second system is the lamel- 
lar cationic lipid-DNA complex structure (Rad ler et at, 19971 
|Lasicefa/., 1997| ( also referred to as lipoplex or phase; 
see below in Fig.0]for a schematic representation) which con- 
sists of a lamellar bilayer stack with one-dimensional arrays 
of DNA intercalated. Experimental evidence suggests that the 
DNA arrays in different galleries are orientationally ordered 
(Battersb y et ai, 19 98 1, implying the presence of DNA-DNA 
interactions across the membrane. Previous theoretical mod- 
eling of this observation was based on the assumption that 
this coupling is mediated through elastic (j Harries et ai, 2003| > 
or electrostatic ( Schiess el and Ajanda-Espmoza^oTH) interac- 
tions. In the latter study, however, it was assumed that the mem- 
brane is infinitely thin, thus neglecting the presence of regions 
with low dielectric constant. Whether direct electrostatic DNA- 
DNA interaction may lead to the coupling is thus still an open 
question. 

The present work analyzes the interaction of a macroion 
with a charged lipid bilayer. Particular focus is put on 
the role of the low dielectric medium inside the mem- 
brane and the electrostatic coupling across the bilayer. The 
level of our analysis will be Gouy-Chapman theory which 
is based on solving the non-linear Poisson-Boltzmann equa- 
tion (Andelman, 1995 1. Another approximation is that be- 
yond the two apposed monolayers of the charged mem- 
brane also the membrane-adsorbed macroion will be mod- 
eled as a very large planar surface. This allows us to ne- 
glect end-effects and, thus, to solve the one-dimensional 
Poisson-Boltzmann equation. While some aspects of elec- 
trostatic interactions across a lipid bilayer have previously 
been analyzed ([Nelson et ai, 19751 Chen and Nelson, 2 000 
|Genet et ai, 2000T " |TanerT^ragTiHmdHi^b05F - mostly on 
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the level of Poisson-Boltzmann theory - we believe our anal- 
ysis is the first that gives analytical expressions for the free 
energy of the large macroion-bilayer system in all interaction 
regimes. 

A second objective of the present work is to ana- 
lyze whether the electrostatic interactions between the DNA 
molecules across the cationic lipid membranes of a lamellar 
lipoplex are energetically strong enough to induce orientational 
ordering of the DNA arrays. Based on our free energy ex- 
pressions for macroion-bilayer interactions, we predict that, 
indeed, electrostatic interactions are expected to significantly 
contribute to the ordering, at least for conditions of low salt 
content (corresponding to a large Debye length). 

Finally, we backup our approximative one-dimensional 
Poisson-Boltzmann model by numerical solutions of the two- 
dimensional Poisson-Boltzmann equation for the lamellar 
lipoplex geometry shown below in Fig. |4] (that is, rigid rods 
intercalated between a perfectly planar membrane stack). Our 
procedure to numerically solve the Poisson-Boltzmann equa- 
tion is a new lattice Boltzmann method which is related to a 
previous application (Horbach and Frenkel, 2001 1. The results 
of our numerical calculations are in close agreement with the 
analytical predictions. 

Electrostatics of a bare membrane 

We start our analysis by considering a charged membrane with 
no macroions present. In this case, it is the two lipid head- 
group regions belonging to the apposed leaflets of the bilayer 
that may interact electrostatically. As outlined in the Introduc- 
tion, the Poisson-Boltzmann model will be used to character- 
ize the electrostatic potential at the headgroups and to calculate 
from that the corresponding free energy. 

To facilitate reading of the subsequent parts, we shall first 
provide a summary of the familiar Poisson-Boltzmann model 
for single charged lipid layer. Formally, this case describes 
an infinitely thick (and thus electrostatically completely decou- 
pled) bare lipid bilayer. 

Single charged lipid layer 

We consider a single, negatively charged, lipid layer. As 
is often the case in experimental systems, the lipid layer is 
mixed, containing a monovalently charged, and an uncharged 
(zwitterionic) lipid species. Denoting the mole fraction of the 
charged species by <j> and assuming (for the sake of simplicity) 
that both lipid species occupy the same cross-sectional area 
a, the surface charge density of the lipid layer is a = —cfte/a 
where e denotes the elementary charge. The minus sign in this 
equation accounts for the negative lipid charges. 

Let us place the lipid layer at the y, z -plane of a Cartesian 
coordinate system such that the headgroup charges are located 
at x = 0. The aqueous region extends along the positive x- 
direction (x > 0) whereas the lipid tails are directed in the 
opposite direction (x < 0) as shown in Fig.^ The subject of 
Poisson-Boltzmann theory is the characterization of the small, 
mobile, ions in the aqueous solution. These ions reflect the 
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Figure 1: Schematic representation of a binary lipid layer, 
consisting of negatively charged and neutral lipids. Some co- 
and counterions are depicted. The graph shows the (reduced) 
electrostatic potential ip(x) of a single charged lipid layer with 
composition <fi = 0.5 and Debye length Id = 100 A. Note that 

into the calculation of ip{ x ) a ls° enter Ib — 7 A, and a = 70 A . 
The surface potential for this case is ^ = —8.28; see Eq.|3] 



presence of salt with bulk concentration n . In thermal equilib- 
rium, a diffuse ionic layer, consisting of co- and counterions, 
is built that partially screens the fixed membrane charges. The 
electrostatic potential <£> can be calculated if the membrane sur- 
face charge density, er, the dielectric constant of water, ew, and 
the salt concentration in the bulk, no, are known. Note that for 
the flat, large, homogeneous lipid layer in Fig. [2 the electro- 
static potential $ = $(x) depends only on the a>direction. 

In Poisson-Boltzmann theory it is convenient to work, in- 
stead of <£>, ew, and no, with related quantities, namely with the 
reduced (dimensionless) electrostatic potential ip = eQ/ksT 
(where ks is Boltzmann's constant, and T the absolute temper- 
ature), with the Bjerrum length Ib = e 2 / (AireoewkBT), and 
with the Debye screening length Irj = 1/k, defined through 
n 2 = 2e 2 no/ '(eoewksT) (where in the latter two definitions 
eo is the permittivity of free space). In terms of these quanti- 
ties, Poisson-Boltzmann theory leads to the Poisson-Boltzmann 
equation ip"(x) = k 2 sinh-0(x) that the reduced potential ip(x) 
has to fulfill subject to the two boundary conditions ip(x — » 
oo) = and ip'(0) = 2np where we have used the short 
notation p = po4> with the constant po = 2ttIbId/o>- (Note 
that a prime denotes the derivative; i. e., ip'(x) = dip(x)/dx 
etc.) Integrating the Poisson-Boltzmann equation once leads 
with %jj{x — > oo) = to 



-2k sinh 



(1) 



After a second integration we obtain the final solution of the 
Poisson-Boltzmann equation for a single monolayer 
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(2) 



e KX coth [arsinh(p)/2] — 1 

which is plotted for one characteristic case in Fig. ^ We 
note that from Eq. [2] one may obtain the surface potential 
\P = xj)(x = 0). The result is 



= —2 arsinh(po<; 



(3) 
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For later use we also note that in terms of the surface potential 
VP the potential ip(x) is given by 



ip(x) = -2 In 



tanh(*/4)\ 
■tanh(*/4) J 



(4) 



Poisson-Boltzmann theory also predicts that the relation be- 
tween the electrostatic potential and surface charge density can 
be used to calculate the free energy per lipid through a charging 
process c Evans and Wennerstrom, 1 994 1 



f m i=a J d~o = -k B T J 4-0) cty (5) 



where <E> = $(x — 0) is the electrostatic surface potential. The 
integration can be carried out because the relation ^(0) is given 
in Eq.|5J the result is 



fml(4>) = 2</> 



P 



ln(p + q) 



(6) 



with q = ^Jp 2 + 1 (recall p = <ppo). Note that here and in the 
following we shall express all energies in units of ksT. Eq.[6] 
requires <f> > 0. Yet, the convention f m i(—(f>) — f m i(4>) en ~ 
sures applicability of Eq.|5]also for positively charged surfaces 
where < 0. Note finally that Eq.|3]implies f' ml {4>) = -*(</>). 

Charged lipid bilayer 

A lipid bilayer consists of two apposed lipid monolayers 
that form a thin film of low dielectric constant. The electro- 
static properties on one side of the bilayer may, in principle, 
depend on those of the other side. This coupling is the subject 
of the present paragraph. 

In the following, we shall refer to the two apposed mono- 
layers of a lipid bilayer as external and internal) and denote 
their compositions by 4>e and </>/, respectively. Fig. |2 illus- 
trates 1 a number of additional structural properties: The lipid 
bilayer has hydrophobic thickness d with the charged head- 
groups of its two monolayers being located at x = and 
x = — d, and the dielectric constant in the inner hydrophobic 
region (— g? < x < 0) is e^. There are three different regions 
in which we denote the reduced electrostatic potential, ip(x), 
by ipi(x) (for — oo < x < —d), by iPl(x) (for — d < x < 0), 
and by iPe{x) (for < x < oo). The reduced potential then 
fulfills the Poisson-Boltzmann equation in the aqueous regions 
adjacent to each monolayer and the Laplace equation inside the 
membrane 



ip'i(x) — n 2 svahipi 

4>'l(x) = o 

i^e{x) — k 2 sinh tpE 



x < —d 
- d < x < 
< x 



It is convenient to define a coupling parameter 
( |Winterhalter and Helfrich, 1992> through 



H 



e L l 



D 



e\vd 




\|/ L (x) 



Figure 2: Schematic illustration the lipid bilayer, located 
within the region —d < x < 0. The reduced (dimensionless) 
electrostatic potentials are ^i(x) and iPe(x) m the aqueous re- 
gions, and iPl(x) between the charged monolayers, each of sur- 
face charge density, <tj — —e<fii/a and <je — —c^e/o,- The 
reduced surface potentials are denoted by tyj and ^>e- 

With that, the boundary conditions at the two interfacial regions 
read 



V4(0) - HdKiP' L (0) = Ikp^e 
-ip'l(-d) + HdKip' L (-d) = 2kpq4>i 



(9) 



Note that, as we expect, for H = the two conditions reduce 
to the corresponding ones for two single, electrostatically de- 
coupled, monolayers. The potential between the two charged 
monolayers fulfills a one-dimensional Laplace equation and is 
thus given by the linear function, 



(10) 



where = ipi(-d) = ip L {-d) and ^ E = ^z(0) = ^s(0) 
are the surface potentials at the internal and external monolayer, 
respectively. We thus have ip' L — A ^ /d where 



(11) 



is the difference between the reduced electrostatic potentials 
across the bilayer. From the first integration of the Poisson- 
Boltzmann equation (see Eq. Q and the boundary conditions 
(see Eqs.|9} we find the surface potentials 



^ e = — 2arsinh(po0_E) ^ i = — 2arsinh(p </>/) 

where we have introduced the effective compositions 



iJA* 
2po 



2 P q 



(12) 



(13) 



(7) We note that the surface potentials, "f^ = ^e{4>e, <t>i) and 
^e — ^e(4>e, <^j), depend on both monolayer compositions, 
<pE and 4>i, because of the electrostatic coupling between the 
charged monolayers. Combination of Eg. 1121 and Ea.ll3lleads 
to A* = fLi(^i) ~ /Me) or, equivalently, 



A^ 

(8) — — = arsinh 



!>j - — A* 
2 



arsinh 



Po<t>E + yAf 
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(14) 
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which represents a transcendental equation for the potential dif- 
ference A\& across the membrane (Baciu and May, 20041 . Of 
course, for a symmetric bilayer, where 4>e = <j>i, the solution 
ofEq.dis A* = 0. 

The electrostatic free energy, fu (measured per lipid pair) 
is given by a charging process analogous to Eq.|5] 



-W 



-W 



*x(0x,O)dfo 



(15) 




where here we first charge the internal monolayer, and then the 
external one (any other path of charging would lead to the same 
result). We obtain 



(16) 



The first two terms represent the charging energies for two sin- 
gle monolayers of effective compositions <\>e and The last 
term, H A^> 2 /4p = e i aA$ 2 /2d is the contribution of the ca- 
pacitor between the charged monolayer interfaces where A$ 
is the potential difference between the external and internal 
monolayer. 

For H — there is no electrostatic coupling between the 
charged monolayers. In this case, fu — f m l(<pE) + fmi(4>i) is 
given by the sum of the individual free energies of two isolated 
monolayers. 



Macroion-bilayer interaction 

After heaving dealt with a bare lipid bilayer, we are now in 
the position to add a positively charged macroion to the exter- 
nal side of the bilayer. To be specific, we place the charges 
of the macroion at position h along the a;-axis as shown in 
Fig- below. It is convenient to express the charge density, 
up = cf)pe/a, of the macroion in terms of the compositional 
variable <f>p (with < 4>p < 1). For 4>p = <f>E the charge 
densities on the macroion and on the external monolayer match 
each other; that is, they would be opposite in sign and equal in 
magnitude. 

As for the isolated bilayer, the potentials, ipi(x), iPl(x), 
and iPe{x), are determined by Eqs. Q The only difference 
is that the last one of Eqs. [7] applies only within the region, 
< x < h, between the external monolayer and the macroion. 
The boundary condition at the macroion is ip' E (h) = 2Kpo<fip. 
The other two boundary conditions (see Eq. [9} remain un- 
changed. The surface potentials = ipi(-d) = i/jl(—cI), 
^ E = 4>l{0) = ^e(O), and * P = ijj E (h) depend on 
the three compositions, <pi, <pE, and (pp. We thus write 

= #r(0i, 4>eAp), = ^e^iAe^p), and * P = 
typ(4>i, (f>E> 4>p)- Fig.|3]schematically illustrates the macroion- 
bilayer system. The free energy of the macroion-membrane 
system, measured per lipid in one monolayer, is again given by 
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Figure 3: Schematic illustration the lipid bilayer, interacting 
with a macroion. The difference to Fig.|2]is the presence of the 
positively charged macroion of surface charge density up = 
ecj>p I a at position x — h. The surface potential at the macroion 
is denoted by "Pp. (Depicted is the case of a weakly charged 
macroion.) 

the familiar charging process analogous to Eq.1151 

4>i 4>e 

*/(0/,o,o)<% -J y E {<i>i,4>E,o)d^E 

ii 



+ / ^p(^i,4>E,4>p)d4>p 
" 



(17) 



It corresponds to charging first the internal monolayer (first 
term in Eq. Mil , then the external monolayer (second term in 
Eq.ll7K and finally the macroion (last term in Eq. I17li. We re- 
peat that the final result for / does of course not depend on the 
order in which the surfaces are charged. 

In the following, the distance h between the macroion 
and the bilayer is chosen such that the free energy / 
is minimal. In this case, the potential iPe(x) has the 
same functional form as for an isolated charged monolayer 
( Parsegia n and Gingell, 1972) : for small <pp the potential is 
that of the isolated lipid bilayer, and for large <f>p the potential 
is determined solely by the macroion charges. The small and 
large <\>p cases are separated by another, intermediate, regime 
where the macroion contacts the bilayer (implying h — 0). In 
the following, we shall consider the three possible scenarios 
separately. 

Weakly charged macroion 

If <pp is sufficiently low the potential ip(x) is equivalent 
to that of an isolated bilayer within the region — oo < x < h. 
The surface potentials at the membrane are thus, again, given 
by Eq.^J and the value of the potential at the macroion, 

* P = -2arsinh(p o 0p) (18) 
is negative. Recalling the functional form of Eq. 0] and 
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ip E (h) — it follows that 



defined by the condition 



tanh(* P /4) \4> E ) \^T+Jp^f -I 



e nh = tanh (* £ /4) 



(19) 

which allows us to calculate the equilibrium distance, h, be- 
tween the macroion and the external monolayer. With growing 
<f>P the distance h decreases until close contact is established 
for a particular <\>p which fulfills the relation <fip — 4>e- The 
condition 

0p < <t>E + t^-A* (20) 

thus defines the weakly charged macroion regime. Recall that 
the potential difference, A\&, fulfills Eq. [21 The presence 
of the oppositely charged macroion reduces the free energy of 
the isolated bilayer by the charging free energy of the isolated 
macroion 



/ = fml{4>E) + f m l{h) + J^-A* 2 - f ml (^ P ) 



(21) 



Highly charged macroion 

For sufficiently high <pp, the potential between the macroion 
and the external monolayer (0 < x < h) is equal to that of the 
isolated macroion, implying 

, , s „, (e- K(x - hs i -tanh(*p/4)\ 

tb E (x) = -21n -. — r> (22) 

l^e-^- ft) + tanh(*p/4) y V 7 

Hence, both the surface potentials at the macroion and at the 
external monolayer are positive, 

= 2arsinh(po^p) = 2 arsinh(po</>p) (23) 

whereas the surface potential of the internal monolayer 



-2arsinh(po</>/) 



(24) 



can be either negative (cf>] > 0) or positive (<f>j < 0), depend- 
ing on the sign of (f>j. Inserting the definitions for (f>j and <f>E 
(see Eqs.ll3> we find a transcendental equation for the potential 
difference A^ across the bilayer 



A* 



= arsinh 



Po> 



bi A* 

2 



arsinh 



Po<pE + yA* 



(25) 

Because of ^>e = V'b(O) according to Eq.[22] and with Eq.1231 
we can determine the equilibrium distance h through 



tanh(*p/4) 



tanh(* B /4) \<j> P J \^/TT(p^e) 



y/1 + (p o 0p) 2 - 1 



1 



(26) 

With decreasing <pp the distance h decreases too until close 
contact is established for a particular <fip which fulfills the rela- 
tion (ftp = cf>E- Hence, the highly charged macroion regime is 



9P > 0E 



H 



A* 



(27) 



where, the potential difference, A 'J, results from the solution 
of Eq.[25] Note that the surface potential at the internal mono- 
layer, fyj, is positive if HA^ > 2p (f>j, and combination with 
Eq.l27lreveals that <f)p > cfip + 4>i is a necessary condition to 
reverse the sign of both surface potentials, / and ^ e- 
The presence of the macroion leads to the free energy 




(28) 



Intermediate case: Close contact 



In between the weakly and highly charged macroion 
regimes (defined by Eq. |20| and Eq. I27> the macroion and the 
external monolayer are at close contact. Here, the surface 
potentials at the lipid bilayer are determined by the equations 



\&j = —2 arsinh \po(4>i 
^P_b = —2 arsinh [po{<pi 



- <Pp)} - 



2po 
H 



(29) 

b E - 4>i) 



Note that because of h = it is *S>e = "Jp- The free energy 
can be calculated by the usual charging process 



, 4>e - 4>p = o) d4>! 




^e{(/>i,4>) d4> 
(30) 

which corresponds to first charging the internal monolayer and 
then the together the external layer and the macroion. We ob- 
tain the result 



(31) 



which is the sum of the electrostatic free energy of a single 
isolated monolayer of composition \<f)j + <f>E — <fip\ and the 
capacitor energy H A^ 2 /4po- Note the particular case <j>p = 
4>e + <fii for which / = po(f)j /H. 

Two limiting cases 

The free energies / in Eqs . 12 1 1 l28l and I3T1 cover all possible 
interaction regimes of a macroion with an oppositely charged 
lipid bilayer. The actual calculation of /, however, requires 
to solve an additional transcendental equation for A 'J. There 
are two limiting cases for which / can be expressed directly in 
terms of the membrane and macroion charge densities. One 
is the Debye-Hiickel limit which is based on linearizing the 
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Poisson-Boltzmann equation. The other one corresponds to a 
very thin (transparent) lipid bilayer for which the limit d — > 
may be considered. These two cases are discussed in the fol- 
lowing. 

Debye-Hiickel limit 



Of course, at the transition from one regime into the other the 
free energy is continuous (but not necessarily smooth). That 
is, for 4>p = 4>e — {4>e — 4>i)H/(l + 2H), the free energy 
/ in Eq. [35] is equal to that in Eq. |40| Similarly, for <\>p = 
4>e + H(4>e + 4>i), the free energy / in Eq.[39]is equal to that 
inEq.gO] 



For sufficiently small charge densities the Poisson-Boltzmann 
equation can be linearized. Formally, the Debye-Hiickel case 
may be viewed as the limit corresponding to Id — ► (or, 
equivalently, po — ► 0). In this case, all three regimes of 
the interacting macroion-bilayer system can be calculated 
explicitly. 

In the weakly charged macroion regime linearization of 
Eq-Elgives rise to 



A* 



fa - <h 
1 + 2H 



(32) 



We can use this result to calculate the condition for the system 
to be in the weakly charged macroion regime (see Ea.l20> 



bp < <pE 



H 



1 + 2H 



r) 



(33) 



as well as the relation for the distance, h, between the macroion 
and the external monolayer (see Eq.ll9> 



„ k h 



H( 



<f>p{l + 2H) 
and the free energy (see Ea.l21> 

Po 1 + 2H 



(34) 



(35) 



Similarly in the highly charged macroion regime, lineariza- 
tion of Eq.[25]results in 



A* = 2 Pq (^e + <t>i) 



(36) 



This result is used to calculate the condition for the system to 
reside in the highly charged macroion regime (see Ea.l27> 



fa>fa + H(fa + fa) 



(37) 



as well as the relation for the distance, h, between the macroion 
and the external monolayer (see Ea. l26t 



„Kh 



fa + H{(f>E + 4>l) 

and the free energy (see Eq.l28> 



]_ 

Po 



b 2 i ~ <Pe 



-HI 



+ fa 2 



(38) 



(39) 



Finally, for the intermediate case of close contact, Eq. 
already provides an explicit expression which in the Debye- 
Hiickel limit reads 



f_ 

Po 



+ ^ [fa 



(40) 



Transparent bilayer 

For very thin membranes (d — * 0) or in the limit of 
large Debye length (Id — > oo) it is H — > oo, and the 
membrane appears transparent to the electric field. 

Consider first the low macroion charge limit. For H — > oo 
we obtain from Eq.^Jthe potential difference 



A* - (fa -<t>i)^0 



(41) 



The effective compositions, defined in Eq.^J are then simply 



bE = 4>i = - {fa + <t>i) 



(42) 



We also see that the condition for residing in the weakly 
charged macroion regime becomes (f>p < <fi where we have 
defined the average composition of the lipid bilayer <\> = (cf>i + 
fa)/ 2. The distance between the macroion and the external 
monolayer is again given by Eq.[^]with fa = <f>. Finally, the 
free energy (see Ea.l21> becomes 



/ = 2f m i((j)) - f m l{fa) 



(43) 



For (f>p > 4> there is close contact between the external mono- 
layer and the macroion. A third, high macroion charge regime, 
does not exist for a transparent membrane. In the close contact 
regime, Eq. ^predicts / = / TOi (0j + <\> E - P ). 

Electrostatic coupling between the DNA 
monolayers in lamellar lipoplexes 

In this section we shall apply our results for the electro- 
static coupling to a specific system, the lamellar lipoplex. As 
mentioned in the Introduction, there is experimental evidence 
for the orientational alignment of the DNA arrays across the 
cationic bilayers (Batt ersby et ah, 1998) . Such an alignment 
implies some kind of energetic coupling. Our results from 
the previous sections enable us to calculate the electrostatic 
contribution to this coupling. Fig. [4] depicts 2 the structure of 
the lamellar lipoplex with its one-dimensional arrays of DNA 
monolayers being intercalated between a stack of cationic bi- 
layers. 

Our objective is to use the formalism developed in the pre- 
vious sections to (roughly) estimate the interaction strength of 
the DNA molecules through a charged lipid bilayer. To do 
so, it is convenient to consider only a single, perfectly planar, 
cationic lipid bilayer with two oppositely charged DNA arrays 
attached to each of its apposed monolayers. Assuming all DNA 

2 fig 4 here 
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Figure 4: Schematic illustration of the lamellar lipoplex struc- 
ture. The (double stranded) DNA molecules, depicted as 
long rods, form one-dimensional arrays between each pair of 
cationic bilayers. Arrays in consecutive membrane galleries are 
orientationally aligned and shifted "out of phase". Some lipids 
are schematically shown; d denotes the membrane thickness. 



strands to be parallel we shift one array with respect to the other 
one. The corresponding free energy is expected to adopt its 
minimum and maximum for the two DNA arrays being "out 
of phase" and "in phase", respectively (see states A and B in 
Fig.0. 

At this point it should be noted that the recent theoreti- 
cal prediction of a so-called sliding columnar phase is based 
not only on an energetic penalty of laterally shifting the two 
DNA arrays (as shown in Fig.[5J but also on a second energy 
contribution associated with angular changes of the DNA ar- 
rays with respect to each other i jO'Hern and Lubensky, 1 998 
IGolubovic and Golubovic, 1998) . Clearly, our present simpli- 
fied model allows to estimate the former but not the latter con- 
tribution. 

We calculate the corresponding free energy difference, AF, 
between states A and B per persistence length, £ = 500 A, of 
double stranded DNA. Let us denote by N the number of lipids 
that interact electrostatically with the DNA. There is some arbi- 
trariness in choosing N because the planar macroion geometry 
(see Fig.|5jl is not a very realistic model for a DNA molecule. 
Still a reasonable estimate is N = 2£i?o/a w 150 where 

Rd = 10 A is the radius of the DNA rod (typically, a « 70 A is 
the cross-sectional area per lipid). The additional factor of two 
in N = 2£Rd/ a results from the two monolayers that interact 
with a single DNA molecule. 

The main difference of states A and B in Fig. [5] is the sym- 
metry of state B with respect to the bilayer's midplane. That 
is, for state B the potential difference across the bilayer A^ff 
vanishes (corresponding to H = 0) and we can write for the 
difference in free energy between state A and B 

AF = F A — Fb = N[f(4, H) - f(<j>, H = 0)] (44) 

where / is the free energy per lipid in one monolayer, and is 
given by Eq.|^for low macroion charge, by Eq.|^for close 
contact, and by Eq.|28]for high macroion charge. 

To specify the compositional variable <j> in Eq.|^]we note 
the ability of the lipid bilayer to adjust its local lipid composi- 
tion. That is, close to the DNA adsorption site the local lipid 



(A) (B) 




D L 



Figure 5: Our structural model to calculate the degree of elec- 
trostatic coupling of two DNA rods (depicted with rectangular 
cross-section) across a cationic bilayer: in state (A) the DNA 
rods are "out of phase" so that potential gradients along the 
membrane normal direction can form. In state (B) the DNA 
rods are "in phase", implying transmembrane symmetry and 
thus vanishing gradients. The difference in free energy of states 
(A) and (B) is a measure for the degree of electrostatic cou- 
pling. Note that, in general, the mixed cationic membrane may 
optimize its composition (mole fraction of charged lipids); <j>D 
and (pL denote the compositions at the DNA adsorption region 
and at the bare membrane, respectively. 



composition, which we shall denote by (f>jj, may differ from 
that at the DNA-free bilayer region (denoted by 4>l in the fol- 
lowing). The ability for local demixing can (approximatively) 
be included in our calculation of AF. In order to derive an 
explicit relation we shall assume small electrostatic coupling 
H <C 1 for which we can expand AF up to first order in H. 
The result is 

AF = -Sa*(&>, M} 2 = ~ \SMd) T fUMf 
4p 4p 

(45) 

where the minus sign refers to the low and the plus sign to the 
high macroion charge regime (in the small H limit the small 
macroion charge regime is <f>p < <pp> and the high macroion 
charge regime is <pp > <\>p,\ the close contact regime reduces to 

4>p = 4>d)- 

Let us analyze the consequences of Eq.|^] For frozen lipid 
composition, where <pn = <f>L = 0, we find AF = in the 
low macroion charge regime. This is expected because for low 
4>p the macroions do not affect the membrane potentials. Yet, 
the membrane is symmetric, implying A^> = 0, and the elec- 
trostatic coupling across the bilayer is irrelevant. 

Consider now the high macroion charge regime for frozen 
lipid composition (<firj = <\>l — 4>)- Here Eq.l45lbecomes 

TT 

AF = -AN— arsinh 2 (p o 0) (46) 
Po 

In Fig. |6] we plot AF according to Eq. |46] (with H = 

CLlo/twd and p — 2-kIdIb/o) for Ib = 7 A, a = 70 A , 
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N = 150, e w /e L = 40, and d = 30 A as a function of the De- 
bye length Id- Curve (a) is derived for = 0.5, curve (b) for 
4> = 0.25. In both cases AF can become considerably larger 
than ksT under low salt conditions. 




20 40 60 80 100 



Id/A 

Figure 6: AF according to Eq.|46]for l B = 7 A, a = 70 A 2 , 
N = 80, e w /e L = 40, and d = 30 A as a function of the 
Debye length Id- The two curves correspond to <f> = 0.5 (a) 
and <fr = 0.25, (b). The broken line shows AF according to 
Eq.g7]for cj) P = 0.5. 

Let us also estimate the influence of local lipid demixing. 
To this end we consider a weakly charged membrane. The 
presence of the DNA's induces accumulation of charged lipids 
near the DNA. From an electrostatic point of view, local charge 
matching {tpD « </>p) would be the optimal situation. The 
matching 4>d = <f>p would be opposed by the demixing free 
energy of the lipids but we may still consider the limiting situa- 
tion that all charged lipids are driven towards the DNA adsorp- 
tion site, implying <\>d = 4>p and </>l = 0. Eg. 145 1 then reads 

AF = —N — arsinh 2 (>o<M (47) 
Po 

Comparison of Eq. [46] with Eq. &\ shows that the lateral mo- 
bility of the lipids tends to reduce the electrostatic coupling 
of the DNA arrays through the cationic bilayer. Sufficiently 
strong electrostatic coupling is generally found only if the De- 
bye length is considerably larger than that under physiological 
conditions where Id ^ 10 A. 

Solving the Poisson-Boltzmann equation 
using the lattice Boltzmann method 

For our numerical calculations we developed a lattice Boltz- 
mann method who's stationary state solves the Poisson- 
Boltzmann equation in the unit cell of the lamellar lipoplex 
geometry. The lattice Boltzmann method simulates densities 
fi on a lattice. Each of these densities is associated with a ve- 
locity v; and they evolve with (scaled) time t through the lattice 
Boltzmann equation 

fi{v + v h t + 1) = fi(r,t) + -L[/°(r,t) - fi( T ,t)]. (48) 



d 




Figure 7: Simulation results for the potential ty(r) in a lamellar 
lipoplex. The two diagrams correspond to the "out of phase" 
(A) and "in phase" (B) cases (see also Fig.|5}. Contour lines 
are drawn at integer values of ^. The simulation is carried out 
for a 1 12 Ax 160 Asystem. The membrane width is d = 30A, 
the DNA radius is Ro = 10A, and the distance between the 
DNA and the membrane is 3A. We use l D = 100A, 
7 A and e^/ew = 1/40. The membrane is homogeneously 
charged with composition <fi = 0.25 of monovalent cationic 
lipids. The charge density on the DNA is twice as large as 
that on the membrane, amounting to a linear charge density of 
— e/1.7A along each DNA rod. The location of the membrane 
charges is indicated by the vertical dashed lines. 
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We identify the potential \P(r, t) = ^ /j(r, t) and choose the 
distribution ff such that it has the moments 

J2f° = 7 ( r ) 0{*W sinKgQ + 2 7 r/ J3P/ (r) ^ ^ 



(49) 



£/°v*v 4 



*(l,t). 



We further choose r(r) = [e(r)/e^ + l]/2. Summing equation 
(I48> we obtain the macroscopic equation 



dfr(r,t) 
9t 



£(£), 



W(r, t) 1 = 6>(rK sinh (*(r, *))+p/(r) 

(50) 

so that the stationary solutions of this simulation solves the 
Poisson-Boltzmann equation for fixed input charge density 
p/(r) and a given field of dielectric constants e(r). The step 
function 0(r) determines the accessibility of the space to salt 
ions. Further details on the method will be given elsewhere. 

We set up a simulation with two lipid bilayers and two DNA 
molecules imposing periodic boundary conditions. In our sim- 
ulations one lattice spacing corresponds to 1 A. We exclude 
salt ions from the membrane regions and from the interior of 
the DNA (this is achieved by the step function 8(r) in Eg. 15 Oil 
and calculate the free energy AF(x re i) as a function of the 
scaled lateral displacement x re i of the DNA arrays in consecu- 
tive galleries of the lipoplex. A snapshot of the potential ^(r) 
calculated by our simulations is shown in Fig. [71 The resulting 
free energy per unit length is then multiplied with the persis- 
tence length of the DNA, £ = 500 A, to give the change in free 
energy with the displacement x re i in units of fc^T. This result 
is shown in Fig.[8]and is in good agreement with the theoretical 
prediction of Eq.|46] More specifically, Fig.[8]predicts a maxi- 
mal change of AF » —9ksT which happens to coincide with 
the corresponding value of curve (b) in Fig.|6]for Id = 100A. 

To sum up, the validity of our approximate analytical 
model for the electrostatic coupling of the two monolayers in a 
lipid membrane is supported by the numerical solutions of the 
Poisson-Boltzmann equation (based on the lattice Boltzmann 
method). Generally, the electrostatic communication across a 
lipid bilayer is quite weak due to the large drop of the dielectric 
constant inside the membrane. However, as we have shown for 
a lamellar lipoplex, it cannot always be neglected. 

Conclusions 

The electrostatic interaction of a macroion with an oppositely 
charged lipid bilayer depends on the coupling parameter H = 
£lId/ {ewd). In the (hypothetical) limit H = the interaction 
reduces to that between the macroion and a single monolayer. 
Here, the apposed monolayer - separated from the macroion by 
a slab of low dielectric constant - can be ignored completely. 
Differently expressed, for H = the two monolayers of a lipid 
bilayer are electrostatically decoupled. 



AF 




Figure 8: Numerical results (derived for the same parameters 
as those in Fig.Q for the change of the free energy AF with 
the relative displacement x re i of the DNA arrays with respect 
to each other. The free energies corresponding to the systems 
shown in Fig.[7Jare indicated by (A) and (B). 



For non-vanishing values of H the behavior of the 
macroion-membrane complex is influenced by the electrostatic 
interaction through the hydrocarbon core of the lipid mem- 
brane. We have studied this aspect on the basis of Gouy- 
Chapman theory, solving a one-dimensional Poison-Boltzmann 
equation. (Thus, all approximations inherent in the Poison- 
Boltzmann formalism are also present in this work). We find 
interesting and non-trivial behavior for non-vanishing H, char- 
acterized by three different regimes: A single weakly charged 
macroion does not affect the surface potentials on either one 
of the two membrane leaflets. As the surface charge of the 
macroion is increased, a regime of close contact will be passed, 
leading ultimately to a third regime of high macroion charge 
where the sign of the electrostatic membrane potential, either 
on one or on both membrane leaflets, is reversed. Within our 
one-dimensional model, all characteristic quantities such as 
membrane potentials, macroion-membrane distance, free en- 
ergies, and boundaries between the different regimes can be 
calculated analytically. 

Typically, the ratio between the dielectric constants inside 
the membrane (e^) and in the water phase (ew) is very small, 
£l/cw ~ 1/40. Thus, only for large ratio between the De- 
bye length (I o) and the membrane thickness d is H not small 
compared to 1. Yet, even for H <C 1 electrostatic interac- 
tions across a lipid membrane can be relevant. To this end, 
we have estimated the electrostatic interaction between the or- 
dered DNA arrays in different galleries of the lamellar 
lipoplex structure. For sufficiently large Debye length Id we 
find a significant electrostatic contribution to the experimen- 
tally observed inter-locking of orientationally ordered DNA ar- 
rays across different galleries. 

Finally, the present work employs the Lattice-Boltzmann 
method to solve the Poisson-Boltzmann equation for the lamel- 
lar lipoplex structure. Our numerical calculations account 
for the structural details of the L„ lipoplex and compare well 
with the analytical results derived on the basis of the one- 
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dimensional Poisson-Boltzmann model. We thus believe that 
the one-dimensional model is a useful tool to estimate interac- 
tions between macroions across a lipid bilayer. 
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